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This work considers a server that processes J classes using the 
generahzed processor sharing disciphne with base weight vector a = 
{ai, . . . ,aj) and redistribution weight vector = (/3i, . . . , The 
invariant manifold M of the so-called fluid limit associated with this 
model is shown to have the form M = {x £ -.Xj =0 for j € S}, 
where S is the set of strictly subcritical classes, which is identified 
explicitly in terms of the vectors a and /? and the long-run aver- 
age work arrival rates 7j of each class j. In addition, under gen- 
eral assumptions, it is shown that when the heavy traffic condition 
'^'j^i 7j = X/^=i '^3 holds, the functional central limit of the scaled 
unfinished work process is a reflected diffusion process that lies in M. 
The reflected diffusion limit is characterized by the so-called extended 
Skorokhod map and may fail to be a semimartingale. This generalizes 
earlier results obtained for the simpler, balanced case where jj = aj 
for j = 1, . . . , J , in which case M = R+ and there is no state-space 
collapse. Standard techniques for obtaining diffusion approximations 
cannot be applied in the unbalanced case due to the particular struc- 
ture of the GPS model. Along the way, this work also establishes a 
comparison principle for solutions to the extended Skorokhod map 
associated with this model, which may be of independent interest. 

1. Introduction. 

1.1. Background and motivation. Generalized Processor Sharing (GPS) 
is a scheduling discipline that is used to share a single processing or trans- 



Received February 2006; revised February 2007. 

^Supported in part by the NSF Grants DMS-04-06191, DMI-03-23668-0000000965, 
DMS-04-05343. 

AMS 2000 subject classifications. Primary 60F05, 60F17; secondary 60K25, 90B22, 
68M20. 

Key words and phrases. Heavy traffic, diffusion approximations, generalized processor 
sharing, fluid limits, invariant manifold, state space collapse, queueing networks, Sko- 
rokhod problem, Skorokhod map, extended Skorokhod problem, nonsemimartingale, com- 
parison principle. 

This is an electronic reprint of the original article published by the 
Institute of Mathematical Statistics in The Annals of Applied Probability, 
2008, Vol. 18, No. 1, 22-58. This reprint differs from the original in pagination 
and typographic detail. 



1 



2 



K. RAMANAN AND M. 1. REIMAN 



mission resource among traffic from several sources. Given a single server 
that can process one unit of work per unit of time, and that is being shared 
byJ(l<J<oo) sources or, equivalently, classes, the information needed to 
implement the GPS policy is contained in the weight vector a = (ai, . . . , aj). 
When all classes have a backlog of work, class j is allotted a fraction aj of 
the total capacity of the server. When some classes achieve no backlog by 
using less than their allotted capacity, the remaining service capacity of the 
server is split among the other classes in proportion to their a^'s. A slight 
generalization of this model was considered in [17] and will be considered in 
this paper as well. A more precise description of that model is given here in 
Section 2. More background on GPS is available in [17] and the references 
cited there. 

In [17] fluid and heavy traffic diffusion limits were obtained for the GPS 
model. The diffusion limits were obtained for the special "balanced case" in 
heavy traffic — that is, under the assumption that for every i € {1, . . . , J}, ai 
is equal to the long-run average arrival rate 7j of work brought in by class 
i customers. In this paper we consider diffusion limits under heavy traffic 
in the possibly "unbalanced case" — in which we only assume the overall 
condition J2i=i'^i = 1 = J2i=i7i^ without imposing any conditions on the 
relation between the arrival rates and weights for each class. The main appeal 
of the unbalanced case is the allowance of some degree of priority: classes 
with ai > 7j can be seen as receiving relatively higher priority than classes 
with ai < 7i . An extreme example of this (allowed in this paper but excluded 
from the heavy traffic limit in [17]) is = 0, where class i is served only when 
some other class has no backlog. The advantage of GPS with < ai < ji is 
that, while class i receives relatively lower priority, it cannot be completely 
starved, in the sense that it always receives service with rate at least Oj. For 
example, as considered in [1], a "next-generation" Internet handling both 
real-time traffic that requires a high quality of service and best-effort traffic 
that has less stringent delay requirements can be modeled as a GPS system 
where the weight for the real-time traffic exceeds its long-run average work 
arrival rate, while the opposite holds for best-effort traffic. 

1.2. Relation to prior work. In addition to the practical motivation given 
above, the unbalanced case is also interesting because, as elaborated below, 
new methods need to be developed for its analysis. There are currently 
two main approaches to establishing diffusion approximations for multiclass 
queueing networks: (i) the continuous mapping approach, which is applicable 
when the so-called Skorokhod map (SM), which maps the netput process to 
the corresponding unfinished work process, is well defined and continuous on 
all cadlag paths (see, e.g., [4, 5, 20] and references therein), and (ii) the gen- 
eral procedure outlined in the papers [2] and [21], which is applicable when 
the directions of constraint (or, equivalently, directions of reflection) satisfy 
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the so-called completely-5 condition. Both these methods lead to diffusion 
approximations that are semimartingales. The present setting does not fall 
under either category. Indeed, neither is the GPS SM well defined for all 
continuous trajectories (as can be inferred from Lemma 2.4(i), Theorem 3.6 
and Theorem 3.8 of [16]), nor do the GPS directions satisfy the completely- 
S condition [as can be deduced from the relation (3.1) of this paper]. In 
addition, the GPS diffusion limit is not in general a semimartingale. 

In the balanced GPS case, this problem was circumvented in [17] by the 
use of the so-called extended Skorokhod map (ESM). The ESM is a general- 
ization of the SM that allows for constraining terms that are of unbounded 
variation and therefore enables the pathwise construction of reflected dif- 
fusions that are not necessarily semimartingales [16]. The unfinished work 
process in the GPS model can be expressed as the image of the correspond- 
ing netput process under the associated GPS ESM (see Theorem 4.3 and 
Lemma 4.4 of [17]). Moreover, as shown in Theorem 3.6 of [16], the GPS 
ESM, in contrast to the GPS SM, is well defined and Lipschitz continuous 
on all cadlag paths. In the balanced case, the fluid limit of the netput pro- 
cess is identically zero and, thus, a continuous mapping approach using the 
ESM, instead of the SM, can be applied to obtain diffusion approximations 
for the balanced GPS model (see the proof of Theorem 4.14 in [17] and also 
Section 5.1). 

The situation in the unbalanced case is considerably more complicated 
because the fluid limit of the netput process is nonzero. As a result, charac- 
terization of the diffusion approximation in this case requires a better un- 
derstanding of the cumulative idleness processes associated with each class. 
A similar generalization was considered in the context of open single-class 
queueing networks in [4]. The analysis in [4] hinged on certain properties 
of the SP associated with the queueing network studied in [4], including an 
explicit decomposition of the constraining term of the SP and continuity of 
the mapping that takes the netput process to the cumulative idleness pro- 
cess. As discussed in Sections 4.2 and 5.1, these properties are not satisfied 
in the GPS model and have consequences for both the fluid and diffusion 
analysis. Thus, new techniques need to be developed for the analysis of the 
unbalanced case. 

1.3. Main results and outline of paper. The first step toward establishing 
a diffusion limit is typically the characterization of the long-time behavior 
of the fluid limit (see, e.g., [2] or [4]). In Section 4.2 we explicitly identify 
the invariant manifold for the fluid limit of the GPS model. The "standard 
definition" of an invariant manifold given, for example, in [2] (see also Def- 
inition 5.2 in [13]) is restricted to subcritical fluid limits. Our definition is 
slightly more general in that it also applies to supercritical fluid limits, in 
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which case M provides information on how the fluid hmit trajectories es- 
cape to infinity (see Remark 4.7 for further discussion of this issue). In the 
subcritical case, our invariant manifold can equivalently be described as the 
set of invariant points for the fluid limit, which coincides with the standard 
definition. Specifically, we show that A4 = {x ^ M."^ : = for i G 5}, where 
5 C {1, . . . , J} is the set of classes for which the long-run average work ar- 
rival rate ji is strictly less than the available long-run average service rate 
(after reallocation of service by the GPS discipline), which we refer to as 
the set of strictly subcritical classes. Although it is intuitively clear that if 
7i < Oj, then i £ S, there may also exist z G 5 for which 7^ > Oi — we provide 
a simple, explicit characterization of the set S. In the balanced case, this 
identification is trivial (there are no strictly subcritical classes), while in the 
unbalanced case, as mentioned above, the identification of S is slightly more 
involved in the GPS model due to the lack of an a priori representation of 
the constraining process in the SP decomposition, in terms of the individ- 
ual cumulative idleness processes. However, we show that, when restricted 
to the fluid limit trajectories, it is possible to identify a certain unique de- 
composition of the constraining term into component processes [see Lemma 
3.1 and Lemma 4.4(i)]. As shown in Theorem 4.6, this turns out to be the 
correct decomposition for the identification of the strictly subcritical classes 
and the invariant manifold A4 (see also Remark 4.5). 

If the convergence of the fluid limit to M is sufficiently fast, one expects 
that M would also precisely characterize the state-space collapse that oc- 
curs in the heavy traffic limit. Indeed, the next main step involves proving 
that the strictly subcritical classes vanish in the diffusion limit, and then 
characterizing the behavior of the remaining, critical, classes. The lack of 
a simple representation for the individual cumulative idleness processes (in 
contrast to the networks considered in [4]) once again complicates this anal- 
ysis. Nevertheless, we overcome this difficulty using two key ideas — namely, 
a comparison principle for the GPS ESM (Theorem 3.2), which may be of 
independent interest, and the introduction of the so-called reduced SP (see 
Theorem 5.6). This culminates in our heavy traffic limit theorem (Theo- 
rem 5.7), which is the main result of this paper. The philosophy behind 
our diffusion approximation is explained in greater detail in Section 5.1. We 
believe that our general approach is likely to be more broadly applicable to 
other models that do not satisfy the so-called completely-5 condition such 
as, for example, networks of stations using the GPS discipline [9] or other 
disciplines that involve a "complete sharing" of service between classes (such 
as, e.g., commonly used work-conserving round-robin disciplines). 

This paper is organized as follows. In the remainder of this section we 
describe our notational conventions. Section 2 contains a detailed description 
of the GPS discipline, a characterization of the unfinished work process, 
and definitions of the SP and ESP, culminating in a representation of the 
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unfinished work process in terms of an SP. The fluid limit of the unbalanced 
GPS model is investigated and a simple characterization of the fluid limit is 
obtained in Section 4. Finally, Section 5 contains the statement and proof of 
the heavy traffic limit theorem for the unbalanced GPS model. The proofs 
of the results in Sections 4 and 5 rely on some basic properties of the GPS 
SP that are presented in Section 3. 

1.4. Notation. We now collect together some of the notational conven- 
tions used in this paper. The set of nonnegative reals, nonnegative integers 
and positive integers are denoted by M+, and N, respectively. Given 
a, 6 E M, a A6 denotes the minimum of a and b and a V 6 denotes the maximum 
of a and b. Vector inequalities are to be interpreted componentwise. The 
standard orthonormal basis in M"^ is represented by {ej, i = 1, . . . , J}, and 
the J-dimensional nonnegative orthant M.'^ is equal to {x G R*^ : x > 0}. Let X 
denote the set {!,..., J}. Given E C M*^, X'([0, oo) : E) represents the space 
of -E-valued right continuous functions with left limits and 'Dc([0, oo);M'^) 
represents the subspace of piecewise constant functions with a finite num- 
ber of jumps. Unless indicated otherwise, we will assume that P([0,cxd) : E) 
and 'Dc([0, oo); M*^) are equipped with the topology of uniform convergence 
on compact sets (frequently abbreviated to u.o.c). For / G 2?([0,oo) :£"), as 
usual, f{t—) = lims|t/(s). For t E [0,oo), \f\{t) denotes the total variation 
of / on [0, t] with respect to the Euclidean norm | • | on R-^. The composition 
of two functions / and g is as usual denoted by fog. The identity function 
i : [0, oo) — > [0, oo) is such that i(t) = t for all t G [0, oo). Given a set ^ C M*^, 
represents the indicator function of the set A, which is defined on R'^ 
and is equal to 1 on j4 and is elsewhere. In addition, co[A] denotes the 
closure of the convex hull of A, cone[A] is the cone generated by A and A° 
is the interior of A. Given a matrix D, we use D' to denote its transpose. 
If X",n G N, and X are processes with sample paths in T>{[0,oo) : E), we 
write X" X to denote weak convergence of the measures induced by X" 
on P([0,oo) :E) to the measure induced by X. 

2. Model description. In this section we provide a detailed description 
of the GPS discipline, introduce our assumptions on the workload arrival 
process, characterize the unfinished work process, and define the Skorokhod 
and extended Skorokhod problems. 

2.1. The GPS discipline. We consider a single server queueing system 
with J customer classes, where 1 < J < oo. Each customer arriving into the 
system brings in a certain amount of work that is measured in terms of 
the amount of time required to process it using the server's total processing 
capacity, which is assumed without loss of generality to be 1. The server 
processes the incoming work using the GPS scheduling discipline, which is 
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described below. The work of class i customers is stored in the class i buffer, 
which is assumed to have infinite capacity. For i £l and t £ [0, oo), Ui{t) is 
defined to be the amount of work of class i that is in the system at time t. 

We now describe the GPS discipline. For E QT, we define af to be 
the fraction of the capacity of the server that is given to class i when the 
set of empty buffers is equal to E. We assume that the processor is work 
conserving, so that J2ifE '^f ~ ^ when E ^X. In this paper we focus on the 
case when the fractions af are determined in the following manner by two 
weight vectors a G [0, l]'^ and /? G (0, 1]*^ that satisfy Ylii^x^i — Pi = 1- 
Given the weight vectors, for E = (i.e., when no queue is empty), we define 
af = ai, and for £' C X, we let 

\ a^+ ( II "i I ' ioT E, 

1 0, otherwise. 

For all E CT, af > ai for i ^ I \ E and J2i^E — = 1- Thus, 

ai represents the minimum guaranteed rate assigned to class i. Any excess 
capacity is split among the remaining classes in proportion to the corre- 
sponding components of the vector p. The vectors a and P will be referred 
to as the basic and redistribution weight vectors, respectively. The condition 
f3i > for each i € X is required to ensure that the processor is work con- 
serving. On the other hand, we allow ai = for some i €l. This represents 
the case when the ith class is of relatively low priority and only receives 
service when one of the high priority classes (with aj > 0) does not require 
all of its assigned capacity. More discussion on the relation between GPS 
and priority is contained in [17] (page 104). 

2.2. Characterization of the unfinished work process. In this section we 
first introduce the primitive cumulative work arrival process associated with 
the GPS model and then present a characterization of the unfinished work 
process. 

We assume that all processes are measurable functions defined on the 
probability space (0,jr,p). Let betheP([0,oo):M:[)-valued process such 
that Hi{t) represents the cumulative work brought into the system by class 
i in the interval [0, t] . We suppose that the cumulative work arrival process 
H and initial conditions satisfy the following properties: 

Assumption 2.1. 

1. U{0)eRi. 

2. H ^ I?([0,oo) :]R:[) is nondecreasing and piecewise constant with H{0) = 
0. 
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3. H has a finite number of jump points in every finite interval, almost 
surely. 

Under this assumption, it was shown in Lemma 2.2 of [17] that the set of 
equations (2.1) below uniquely characterizes the set of processes (C/, Ie, E C 
T), where, for E (11, lE{t) denotes the amount of time in [0,t] that the set 
of empty buffers is equal to E. For i £l, 

Ui{t) = UiiO) + Hiit)- ^flE{t) and 

Eel: iiE 

(2.1) lE{t) = j ^{E(s)=E}ds forE'CJ, where 

E{s) = {ieT:Ui{s) = G]. 
The busy time process Tj defined by 

(2.2) T,= ^f^E, 

EcX-.i^E 

represents the cumulative amount of service given to class i. Note that in 
the definition and characterization of the unfinished work process, no as- 
sumption is made as to how the service allocated to a class is divided among 
the customers present in that class. 

The ultimate assumptions that we require on the primitive processes are 
quite weak: H must satisfy a functional strong law of large numbers (see As- 
sumption 4.1) and a functional central limit theorem (see Assumption 5.1). 
The simplest concrete example where this happens is when the J component 
processes of H form compound renewal processes (see, e.g., [19], Lemma 2). 

2.3. Definition of the Skorokhod and extended Skorokhod problems. The 
GPS model for a = (3 with stochastic fiuid inputs was analyzed in [7]. It was 
shown there that the mapping taking the so-called netput process to the 
unfinished work can be represented in terms of a Skorokhod problem. Below, 
in Lemma 2.5, we recall the similar representation that was derived in [17] 
for the unfinished work process U associated with the slight generalization 
of the GPS model described in Section 2.1. First, we need to recall the 
definitions of the Skorokhod and extended Skorokhod problems associated 
with the GPS model. 

Roughly speaking, given the closure G of a domain in R-^, a set of allowable 
directions of constraint d{x) associated with each point x G dG and a path 
tjj, the solution to the associated Skorokhod problem defines a constrained 
version 99 of i/' that is restricted to lie in G by a constraint mechanism ^ — ip 
that is of bounded variation and acts in the direction of one of the vectors in 
d{<p{s)) using the "least effort" required to keep (p in G. The solution to the 



8 



K. RAMANAN AND M. 1. REIMAN 



extended Skorokhod problem is a generalization of the Skorokhod problem 
introduced in [16], which relaxes the bounded variation requirement on the 
constraint mechanism. 

The domain of the GPS Skorokhod and extended Skorokhod problems 
is G = M.'^ and the directions of constraint (sometimes also referred to as 
directions of reflection) are characterized in the following manner by the 
redistribution weight vector f3 E (0, l)"^, which satisfies J2iei f^i ~ ^- 
d-J+i = YuiiiX^il ^ and define 

(2.3) d, = e,- P\ ioTi^l. 

The set of allowable directions of constraint at any point x on the boundary 
dG is then given by 

(2.4) d{x) = \ ^ ajdj : > for i E /(x) L 
where, for x E Mi, we define 



(2.5) I{x) 



{iEJ:Xj = 0}, if ^Xj>0, 
J el 

XU{J + 1}, if ^xj=0. 

The set of directions of constraint on the boundary dG of the domain de- 
scribes how service is reallocated when one or more classes have no backlog. 
Since this reallocation is determined solely by the redistribution weight vec- 
tor /3, the description of the GPS Skorokhod problem only depends on (3 
(and not on the basic weight vector a). Thus, even though our model is 
more general than the one considered in [7, 8, 16], the results from [8, 16] 
can still be applied here. For more intuition on the relation between the di- 
rections of constraint and the reallocation mechanism of the GPS discipline, 
the reader is referred to [7, 17]. 

We now provide the rigorous definition of the GPS Skorokhod problem. 
Let (^{x) = d{x) n {x E M"^ : \x\ = 1} and recall that for rj E ^([0, oo) : M"^), 
|7/|(T) denotes the total variation of r/ on [0, T] with respect to the Euclidean 
norm on M"^. 

Definition 2.2 (Skorokhod problem). Let ip e'D{[0,oo) :R-^) withV'(O) E 
M^^ be given. Then (99, rj) solves the GPS Skorokhod problem (SP) for ip if 
ip{0) = V'(O), and if for all t E [0,oo), the following five properties hold: 

1. V9(t)=^(t)+7y(t). 

2. ip{t)£R-l. 
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3. \r]\{t) < oo. 

4- \v\{i)=I[o,t] -^M.)eORi}^l^l(s)- 

5. There exists a measurable 7: [0, 00) — > M"^ such that j{t) G d^(tp(t)) {d\r]\- 



Note that (p is constrained to remain within M;;^, and that t] changes only 
when If is on the boundary dM.'^, in which case the change points in one 
of the directions in d{ip). If {f,ip — solve the SP for ip, then we denote 
(p = and refer to T as the GPS Skorokhod map (SM). 

The values of ^ e ^^([0, 00) : R-^) for which there exists (p £ T>{[0, 00) : M:[) 
such that (p = r(^), is called the domain dom(r) of the SM T. Theorems 
1.3, 3.6 and 3.8 of [16], together, show that the domain dom(r) is a strict 
subset of ^([0, 00) : M"^) that does not include certain paths of unbounded 
variation. Since diffusion paths are almost surely of unbounded variation, 
the SM is thus inadequate for constructing reflected diffusions associated 
with the GPS model. This necessitates the introduction of the so-called 
extended Skorokhod map, first introduced in [16]. Recall that, for A CM"^, 
co[A] represents the closure of the convex hull of the set A. 

Definition 2.3 {Extended Skorokhod problem). Let ip £ P([0, 00) : M"') 
with ^{0) E Ri be given. Then {ip, rj) solves the GPS extended Skorokhod 



Theorem 3.6 of [16] shows that there exists a unique solution to the GPS 
ESP for ah ip G V{[0, 00) : R-^). Analogous to the GPS SM, if {ip, p-ip) solve 
the GPS ESP for ^p, then we denote p = r(^), and refer to F as the GPS 
extended Skorokhod map (ESM). 

Remark 2.4. Lipschitz continuity of the GPS SM and GPS ESM (with 
respect to the u.o.c. topology) on D([0, 00) : M'^) was established in Theo- 
rem 3.8 of [8] and Theorem 3.6 of [16], respectively. In particular, these 



almost everywhere), and 





(2.6) r?(t)-r?(s)Gco |J d{p{u)) . 




ii{t)-v{t-)eco[d{pm. 
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results imply that solutions to the GPS SM and ESM are unique. Moreover, 
Theorem 1.3 of [16] shows that T{ip) = T{'ip) for every ip S dom(r). 

In order to state the Skorokhod representation for the unfinished work 
process, define 

(2.7) Xi{t) = UiiO)+Hi{t)-a^t 
and 

(2.8) Yi{t)=a^t-Ti{t), 
where, as in (2.2), recall that 

Ti{t)= otflEit). 

Note that, by (2.1), we have \J = X^Y . 
The following result is Lemma 3.4 of [17]. 

Lemma 2.5. Let X he as defined in (2.7) and let V he the GPS SM 
associated with the weight vector /? G (0,1]"^. If {U,Ie,E ^T) satisfy (2.1), 
thenU = T{X). 

3. Some properties of the GPS Skorokhod map. In this section we col- 
lect properties of the GPS directions of constraint that are used for both 
determining the long-time behavior of the fluid limit in Section 4, as well 
as for identifying the diffusion limit in Section 5. A characterization of the 
geometry of the GPS directions of constraint is provided in Section 3.1, and 
a comparison principle is established in Section 3.2. 

3.1. Geometry of the directions of constraint. The main reason that the 
GPS SM is not defined on all continuous functions is because the directions 
of constraint di,i € X, are not linearly independent (see Lemma 5.3 for one 
important ramification of this property and [16] for further discussion of this 
issue). Indeed, the directions of constraint satisfy 

J 

(3.1) ^/?fc(l-/?fc)4 = 0, 

k=l 

which is easily verified by direct substitution. Nevertheless, as shown in Lem- 
ma 3.1, the GPS directions of constraint do exhibit a reasonable degree of 
regularity. This lemma is used to enable a simple description of the action of 
the SM on affine trajectories in Lemma 4.4, and is also used in Theorem 5.6 
to establish a certain reduced representation for the GPS SM. 
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Lemma 3.1 (Geometry of the GPS directions of constraint). For ev- 
ery j G I, the vectors {di, . . . ,dj} \ {dj} span the hyperplane H = {x G 
M'' : {x,dj^i) = 0}. Moreover, the following two properties hold: 

(i) H = U/=i Cj , where Cj = ^^^dk : 0k > 0} for j G I. 

(ii) Given any w ^ H, there exists a unique vector 6 G W'^, such that 9j = 
for at least one j £l and 

J 

(3.2) w = Y,0jdj- 

i=i 

In addition, w admits the representation 

(3.3) Wj = - /3jCT 

for j G I, where, for any set E such that {A; G T : 0^ > 0} dl, 

(3.4) {)<a= ^ - — jr=2^-. — ir = -, — ^ IT- 

Proof. The first statement of the lemma fohows from the proof of 
Lemma 3.1 of [8]. In turn, this statement imphes that U/=i C H and that, 
given any w & H, there exist 6*^ G M, A; G {2, . . . , J}, such that vu = J2k=2 ^kdk- 
li 6k>0 for every k £ {2, . . . ,J}, then it automatically follows that w £ Ci. 
Otherwise, choose / G {2, . . . , J} such that 

, . 0k 

I = argmm 



fce{2,...,j} /3fc(l - f3k) ' 



define 9i = and, using relation (3.1), eliminate di in the representation of 
vu to obtain 

W= > \0k-0l-rTZ TT^ \dk. 



k=l,kjtl 



A(i-A) 



In this case, by the choice of I, 0i <0 and all the coefficients of d^jfe ^ I, 
in the above decomposition are nonnegative. This shows that w £ Ci. So, 
in both cases, we have shown that w G U/=iC'j- Since this holds for every 
w £ H, we conclude that H C [Jj^^Cj, which, when combined with the 
reverse inclusion, establishes property (i). 

For property (ii) , fix G and first observe that the existence of ^ G M.'^ 
satisfying (3.2) follows immediately from property (i). Now suppose that 
there exist i,j £l, £ M;]^ with 9i = and 6i G ]R:( with 6j = such that 

J J 
(3.5) w= ^ 0kdk= XI ^kdk- 

k=l,k^i k=l,kj^j 
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If i = j, then the hnear independence of the vectors dk,k G T \ {i} shows 
that 9 = 6. Now suppose that i ^ j. Then, since relation (3.1) imphes that 

substituting this equahty into the first representation for w in (3.5) yields 
J J 



k=l,k^i k=l,k^{i,j} 
J 

E 



k=iM{id} 



f3k{l-Pk) 



, . A(i-A) , 



Comparing the last display with the second representation for w in (3.5), 
the linear independence of the vectors dk,k £l\ {j}, implies that 



/3fc(l-/3fc) 



&k 



' for k 



The fact that 9i > 0, Oj > and /?i(l - - (3j)] > together imply 

that 6j = 0, 6*4 = and, hence, that 9k = 9k for k £ This establishes 

uniqueness of the vector 9 satisfying (3.2). 

To show the last assertion, using the definition of the directions of con- 
straint and the first property of (ii), we obtain 

for j £ I, which proves (3.3) with a as defined in the first equality of (3.4). 
The second equality follows trivially since 9k = for k £ E \ {k : 9k > 0}. 
Summing the last display over all j such that 9j > 0, we obtain 

j:ej>0 \ j.0j>O / \fe:0fc>O ''^k/ 

which proves the third equality in (3.4) when E = {j-9j > 0}. Since (3.3) 
shows that for every j with 9j = 0, wj = crPj, this now shows that the third 
equality in (3.4) holds for any E CT that contains {j ■.9j>0}. □ 
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3.2. A comparison principle. The main result of this section is a com- 
parison principle for the GPS ESM, which may be of more general inter- 
est. Analogous comparison results for the one-dimensional SM on [0, oo) 
and the two-sided SM on [0,a] can be found in [22] and [14], respectively. 
The statement of the comparison principle involves the projection opera- 
tor vr : M'^ — > M.'^ associated with the GPS SP, which is characterized by the 
property that 

Tr{x)=x if a; G Mi, 

(3.6) 

7r{x)£dRi and tt{x) - x £ d{TT{x)) ifx^M;^. 

The existence, uniqueness and Lipschitz continuity of the GPS projection 
operator was established in Theorem 3.8 of [8]. The comparison principle of 
Theorem 3.2 is used in Theorem 4.6 to provide uniform bounds (with respect 
to initial conditions in a compact set) on the time at which the fluid limit 
of the GPS model reaches the invariant manifold. It is also used in Section 
5.2 to establish "state space collapse" for the diffusion limit (see Theorem 
5.6). 

Theorem 3.2 (A comparison principle). Let vr be the GPS projection 
operator and T be the GPS ESM. If D < u, then 

(3.7) 7r(P) < 7r(i/). 

Moreover, if for £ I?([0, oo); M"^), there exists a coordinate-wise non- 
decreasing function A with A(0) > such that ip{t) =^p{t) + A{t) for every 
t £ [0, oo), then 

(3.8) ^Wit) < r(^) (t) for every t G [0, oo) . 

In the following lemma, we first show how (3.8) can be deduced from (3.7) 
for a general class of ESPs. 

Lemma 3.3. Let tt be the projection operator associated with an ESP 
that has a uniformly continuous ESM defined on V{[0,oo);'R'^). Suppose 
that v <v implies 7r(z>) < 7r(i/). Then, for ip,ip £ !D([0, oo); M"^) such that 
there exists a coordinate-wise nondecreasing function A with A(0) > and 
^{t) = tjj{t) + A(t) for every t G [0, oo), (3.8) holds. Moreover, if the SM is 
uniformly continuous and well defined on 2?([0, oo); M"^), then (3.8) holds 
with the ESM replaced by the SM. 

Proof. When G ^c([0,oo);M-^) (recall that Pc([0, oo); R-') is the 
subspace of P([0, oo); M"^) that has piecewise constant trajectories with a 
finite number of jumps), this can be proved using induction. Indeed, if = 
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to < ii < ^2 < tn are the union of the jump points of ip and -0, it follows that 
T{^p){0) = 7r(^(0)),r(V;)(0) = 7r(^(0)) and for /c = 1, . . . , n - 1, 

r(v)(tfc+i) = vr(r(v)(tfe) + ^(tfc+i) - V(ifc)), 

and likewise for T{ip){tk+i) (see, e.g., (28) of [6] for this construction for the 
SM r — it is easy to see that the same construction also holds for the ESM 
r). Now suppose r{'ip){tk) < T{7p){tk)- By the assumption on A, we have 
'^(tk+i) -tpitk) < Tpitk+i) -i^itk), and thus, (37) implies _that T{'4}){tk+i) < 
r (-(/>) (tfc+i). Since A(0) > 0, (3.7) ensures that T{'ijj){tQ) < r(^)(to) and thus, 
by induction, (3.8) holds whenever & Pc([0, oo); M"^). 

In order to extend the result to general Tp,Tp,A in P([0,oo);]R'^) that 
satisfy the assumptions of the lemma, we need only use the fact that ip, tpjA 
can each be approximated in the uniform norm by corresponding sequences 
and {A('^)} of piecewise constant trajectories with a finite 
number of jumps such that -0^"^ = "ip^^^ + A("\ where A(")(0) = and A*^") 
is coordinate-wise nondecreasing. Indeed, consider the sequence ^„ = {0 = 
to < ti < • • • < f^^, kn £ N} of partitions of [0, oo), such that the nth partition 
contains all points at which ip or tp' has a jump of magnitude greater than 
1/n and the mesh size maxj=i ^.^ |t" — — > as n — > oo. Define 

/(")(t) = /(t,_i) fortG[d,tr), i = l,...,K, 

for / G {^p,tp,A}. Then for every n G N, /^"^ — > / u.o.c. for / G {ipjtpjA} 
and V^"^ = V'^"^ + A(") with A(")(0) = 0, A^") > and A^") coordinate-wise 
nondecreasing. Thus, the inequality (3.8) is satisfied with ip, ip replaced by 
'(/'("^ and 'ip^'^\ respectively, for every n G N. Taking limits as n — > oo, the 
uniform continuity of T then ensures that the inequality (3.8) also holds for 
ip,''P- If the SM is defined and uniformly continuous on P([0, oo); M''), then 
the last statement of the lemma follows due to exactly the same argument 
as that used for the ESM. This completes the proof of the lemma. □ 

Before presenting the proof of the theorem, it will be convenient to in- 
troduce the following notation. Recall that the one-dimensional Skorokhod 
mapping, Fi : ^^([O, oo) : M) "^([0, oo) : M), is given explicitly by 

(3.9) ri(/)(t) = /(t)+ sup [-/(s)]VO. 

se[o,t] 

Also, given j G T, define P^^'^ to be the matrix whose ith column is di for 
i ^ j, and whose jth column is the unit vector ej. Let Q^^^ = I — P^^\ 

Then it follows immediately from the definition (2.3) of di that Q^j^ > for 
hkel, E^eIQ^ = and for every k G T \ {j}, E^eIQi^ = 1- We now 
claim that {Q^-''^)' is the transition matrix of a transient J-state sub-Markov 
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chain, where the chain "dies" after entering state j. To see that it is tran- 
sient, note that for i /i, (Q^-'^)ij = — Pi) > (since, by assumption, 
minjgi/3j > 0). Thus, (Q^^^)' is strictly substochastic and, hence, has spec- 
tral radius a{Q^^^) < 1. We now describe the SP associated with the matrix 
P^^\ Roughly speaking, it has domain M.'^ and direction of constraint in the 
relative interior of the ith boundary face, {x G M.'^. :xj = 0}, given by the 
ith column of P^^^ . More precisely, for j G I, define the set- valued function 
d(j)(-) on dRi as follows: 

{d{x), if Xj > 0, 

\ ^j^j + : aj > for i G J{x) \ , if xj = 0, 

where, for xeRi, J{x) = {i G T : = 0}. Also, let ^(^^^(x) = d^^\x) n{ve 
R"^ : \v\ = 1}. The SP associated with P^^^ is defined as in Definition 2.2, but 
with d^{-) replaced by d'-^^ {■). Let P^^) be the corresponding SM, and let 
vr^-'^ be the corresponding projection operator, characterized by the relations 
(3.6), but with d{-) replaced by d^^\-). SPs of this kind were introduced in 
[11], the results of which show that the properties of P^^^ described above 
guarantee that F^-^^ and tt^^^ are well defined on P([0, oo); M"^) and R*^, re- 
spectively (see also the discussion in Section 2.3 of [8]). 
We now present the proof of Theorem 3.2. 

Proof of Theorem 3.2. In order to prove the first property, fix i>, G 
M"^ such that v <v and define R = 'K{i') and k = '/r(i/). Then the definitions 
of vr and the GPS directions of constraint show that 

(3.10) Yi>i<^ =^ -z>Gd(0) =^ 7r(P) = 0. 

Since tt{i') > 0, this proves the result when J^iei'^i — 0- shall consider 
two cases with J2iei'^i > 0- 

Case 1. J2iei'^i > there exists j £l such that kj > and kj > 0. 

In this case, we first claim that 7r(i/) = ir^^\iy). Indeed, this is a direct 
consequence of the definition and uniqueness of the projection operators vr 
and vr'^-'^ and the fact that kj > implies d^^\ti) = d{K). Since Kj > 0, the 
same argument also shows that 7r(t') = From the definition of the SP 

and the projection operator, it is easy see that T^^\vi) = tt^^\u)l. Moreover, 
the comparison principle proved in Theorem 4.1 of [18] guarantees that 
T^^\i>L){t) < r^^\uL){t) for t G [0,oo), and so, substituting t = 1, we obtain 
the inequality Tr^^\D) < tt^^^v). Thus, TT{i') = ■k^'-'\d) < ir^^^i/) = vr(z^), and 
so the theorem is true in this case. 

Case 2. J2i£i ^ and { j G X : Kj = 0} = {j G T : kj > 0} . 
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We shall argue by contradiction to show that this case cannot occur when 
u < u. For u G M'', let £{1^) = {j ^I:Kj = 0}. Note that the first condition of 
Case 2 implies that —D ^ d{Q) and, hence, k / or, equivalently, £{i>) ^T. 
Next, suppose <?(i^) =X. Then we have k = and so, by the definition of 

^ieT^i — ^- When combined with the ordering v this contradicts 
the first assumption of the case. Now consider the remaining possibility 
when ^{i') CT and <f (t') CX. Then, by the definition of vr and Lemma 3.1, 
w = K — u G d[K) = cone[dj,j G £{i')] Q H. This means, in particular, that if 
9 e M-^+^ is the unique vector in (3.2), then 6'j+i = and {j : Oj > 0} C £{i^). 
Thus, using (3.2)-(3.4) and the fact that i^j = — Wj for j G £{1^), we observe 
that 

(3.11) ^3 = - E ^j = -(i- E /?i)'^<o. 

Since £{i>) / X, an analogous argument also shows that J2jG£{u) ^ 0- More- 
over (3.11), along with the fact that u <v, implies that 

E '^.■< E ^^-^o- 

jee{v) je£{u) 

On the other hand, since £'^{i') = £{i') due to the second condition of Case 
2, we also have 

E E ^.•<o- 

Together, the last two relations imply that Ylj^i^j ^ 0; which contradicts 
the first assumption of the case. This completes the proof of (3.7). 

Since the GPS ESM is Lipschitz continuous, the second assertion of the 
theorem follows from the first due to Lemma 3.3. □ 

4. Long-time behavior of the fluid limit. In this section we consider a 
sequence of GPS systems with an associated sequence {H^} of cumulative 
work arrival processes defined on {Q,,T,P) that satisfy Assumption 2.1. 
Let [/"" and T" be the associated unfinished work and busy time processes 
uniquely characterized by equations (2.1) and (2.2). We also consider the 
associated sequences {-^"} and {5^"}, where and are defined by the 
relations (2.7) and (2.8), respectively, with C/j(0), Hi and Tj replaced by 
J7"(0), Hf and T", respectively. Assume !F is complete with respect to P 
and for n G N, let {Tt} be complete filtrations such that is adapted to 

In Section 4.1 we state the characterization of the fiuid limit of the se- 
quence of unfinished work processes that was obtained in [17]. In Section 
4.2 we identify the invariant manifold for the fluid limit — this constitutes 
the first step toward establishing the heavy-traffic diffusion approximation 
in Section 5. 
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4.1. Characterization of the fluid limit of the unfinished work process. 
Given a sequence {/"} C 'D([0, oo) : K"^), we define the associated fluid scaled 
sequence {/"} C P([0, oo) : R-^) by 

(4.1) 7"(t) = Z!M for tG [0,cx)). 

n 

From Definition 2.2 and the fact that the GPS SM T is Lipschitz contin- 
uous on its domain (see Remark 2.4), it is easy to verify that T is nonantic- 
ipatory in the sense that r{X){t) depends only on {X{s),s < t}. Since ff", 
and therefore X"', is adapted to the filtration {.F"} and since C/" is right 
continuous and f/"" = T{X^) by Lemma 3.4 of [17], this implies that 
is progressively measurable with respect to the filtration {J^"} (see Propo- 
sition 1.13 of [12]). We now assume that the primitive processes satisfy a 
functional strong law of large numbers. Recall that the abbreviation u.o.c. 
represents uniform convergence on compact time intervals. 

Assumption 4.1. 

1. There exists u G M.'^ such that a.s. 

f/"(0) _ 
iim = u. 

n—>oa 11 

2. For each n G N, there exists 7" G M;]^ such that a.s. 

hm ^ = 7"-*, 

m — >oo fji 

where the convergence is u.o.c. 

3. There exists 7 G M;;^ such that 

lim 7*^ = 7. 

n— ►oo 

Recall that — > is the identity map. Define 

(4.2) 2^ = 7-0, 

and note that for j G I, —Vj represents the amount of nominal service ca- 
pacity allocated to class j that is in excess of its mean arrival rate. Also, 
let 

(4.3) X = u + iyi, U = T(X), Y = U -X 
and 

(4.4) T = ai-Y = -ii + u-U. 
The following result was established in [17]. 
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Theorem 4.2 (Fluid limits for the GPS model). Suppose Assumptions 
2.1 and 4.1 hold, and let U and Y be given by (4-3). Then P a.s., as oo, 
U U, Y Y and T T u.o.c. Moreover, if J2i=i 7i = 1 ^'^'^ n = 0, 
then U = and T = ^l. 

Proof. The first statement is Theorem 4.3 of [17]. The condition Ylii=i li — 
1 implies that Ylii=i = ^ and, thus, the second statement follows from 
Lemma 4.4(2) of [17]. □ 

Remark 4.3. In Theorem 4.2 the fluid limit of the unfinished work is 
represented in terms of the GPS SP: U = T{X). As shown in [7], it is also 
possible to equivalently represent the fluid limit as the unique solution to a 
system of coupled differential equations. While the latter may in some sense 
provide a more intuitive description of the fluid limit, it does not extend to 
more general continuous inputs, as is required for the heavy traffic analysis. 
In contrast, the use of the GPS SP (ESP) provides a unified framework in 
which to study the pre-limit, fluid limit and diffusion limit. It is therefore 
natural and more convenient to work throughout with the SP formulation. 
Indeed, this level of abstraction allows one to better understand the con- 
nection between the nature of reallocation of service (as embodied in the 
directions of constraint) and the continuity and monotonicity properties of 
the map that are used in establishing the limit theorems. As a result, we 
expect that this approach may be more readily generalizable to other situa- 
tions, including those in which a simple differential equation characterization 
of the fluid limit is not available. 

4.2. The invariant manifold of the fluid limit. The goal of this section is 
to identify the so-called invariant manifold of the fluid limit. We first con- 
sider the task of identifying the set of strictly subcritical classes, namely, 
the sources whose long-run average arrival rate 7j is strictly less than the 
long-run average service rate available to them (after redistribution of ser- 
vice by the GPS discipline) for all sufficiently large t. For SPs associated 
with other queueing networks that have been studied in the literature (see, 
e.g., [4] or [21]), the linear independence of the associated constraint direc- 
tions helps simplify this task. For example, in the study of open single-class 
queueing networks in [4], the associated unfinished work U* is represented 
as the image of the corresponding netput process X* under the associated 
SM r*: U* = T*{X*) = X* + Y*, where Y* now admits the decomposition 
Y* = Z?^*, where D = {I — P') is a nonsingular matrix, / is the identity 
matrix, P is the so-called routing matrix of the network, the ith column of 
D represents the direction of constraint associated with the face {xi = 0}, 
and ^* is the component-wise nondecreasing process that characterizes the 
cumulative idle time or excess capacity at each station (note that, somewhat 
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unfortunately, the term ^* here is denoted by the letter Y in [4]). The fluid 
limit ^* of the vector of cumulative idleness admits the explicit expression 

= D~^Y , where Y = U — X , with U = T*{X ) being a continuous 
functional of X that is identically zero under the overall heavy traffic con- 
dition (the bar quantities here all refer to the fluid limits of the original 
quantities). In this case, the set of strictly subcritical classes is precisely the 
set of classes j for which 9* = d^*/dt is strictly positive for all sufficiently 
large t. Since 9j can be explicitly recovered from X* , which is itself known 
explicitly in terms of the primitives, this simplifies the determination of the 
strictly subcritical classes. 

In contrast, as shown in (3.1), the GPS directions of constraint are linearly 
dependent and, thus, such a simple linear algebraic relation between the 
netput and the cumulative idle time processes no longer holds in general. 
Nevertheless, using the geometric properties established in Lemma 3.1, we 
show below in Lemma 4.4 that an analogous decomposition into component 
processes is possible for fluid trajectories. However, it is important to note 
that this decomposition does not hold for arbitrary trajectories (see Lemma 
5.3), thus necessitating a more careful analysis of the diffusion limit (see 
Section 5.1). 

Lemma 4.4 (A representation lemma). Given u G and v G M-^, there 
exists e > and x £ I^"^ such that 

T{u + vL){t) = u + xt fort£[0,e). 

Moreover, the following two properties are satisfied: 



(i) There exists a unique vector 9 G M'^'^^ that satisfies 



J+i 

(4.5) X = ^+J2^kdk 

k=l 

and 

(4.6) I^{k€l:9k>0}<^{k€l:Uk = Xk = 0}- 

Moreover, ^j+i > if and only if u = x = and J2j=i < 0- 
(ii) lfu = 0, then x = '?r(z^), e can be chosen to be oo and [J2j=iXj] 

r 



E/=i^,]vo. 



(hi) //E/=iZ^i>0, th 



en 



(4.7) 0< ^ - 



h _ ~ I]fc:6>fc>0 

h - Y.k:e,.>0 Pk 
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and, for every j £ I, 

if e, > 0, 

if = 0. 

Proof. Fix i/ G M-^ and let £ = {k £ I -.Uk = Xk = ^} ■ The existence of 
X G IR"^ and e > results from the fact that images of affine trajectories 
under the GPS SM are piecewise affine (which was proved in Lemma 4.4 
of [17]). Next, note that for all t £ (0,e), the definition of the GPS SP [in 
particular, the relation (2.4) and properties 4 and 5 of Definition 2.2] implies 
that X is the unique vector that satisfies 

{cone[dj,j G £], 
cone[dj,j G2,dj+i], 
\i£=X. 

This immediately ensures the existence of a vector Q G W]^^ that satisfies 

(4.5) and (4.6). To show uniqueness of 0, first taking inner products of (4.5) 
with (ij+i shows that 

(4.10) (x, dj+i) = (y, dj+i) + ej+i. 

If £ ^I, then (4.9) and Lemma 3.1 show that x ~ ^ ^ H , and so ^j+i = 
by (4.10). On the other hand, if £ =1, then (4.10) uniquely determines ^j+i, 
> if and only if J2j£i < and x~''^ ~ (^J+idj+i £ H. The last two 
statements, when combined with Lemma 3.1(ii), establish uniqueness of the 
representation (4.5) and the condition on thus proving property (i). 

Now suppose that = 0. Then a simple consequence of the definition 

(3.6) of the GPS projection is that tt(i')l solves the SP for i^l. Therefore, 
by uniqueness of solutions to the GPS SP, x = '^(^) £ ^+ and e can be 
chosen to be oo. The second relation in property (ii) can be deduced in 
a straightforward manner from the fact that x = ''^i^) ^ and the last 
statement of property (i). The latter also shows that when X)/=i '^j — 0, 
6j+i = and w = x~ ^ ^ H . When combined with (3.3) and (3.4) of Lemma 
3.1, elementary algebra immediately yields property (iii). □ 

Recall from (4.2) that = 7 — a. When X^jgj^j < 0, it is intuitively 
clear that all classes will be strictly subcritical. The interesting case is thus 
when Y^j^x^j — 0- Suppose, in addition, that n = 0, and note that by 
Lemma 4.4(ii), in this case x = ^('^) ^ - Let 9 £ M'^^ be the unique vector 



(4.8) 



Xj 



P3 



+ (3, 



1 Y.k:ek>0 
I]fc:6»fe>0 ^k 



f^k 



1-E 



fc:6»t>0 



Pk 
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in the representation (4.5) for x = '^('^)- With reference to the discussion 
prior to Lemma 4.4, it is natural to introduce the following definition: 

(4.11) 5o(i^)={jGJ:%>0}. 



Remark 4.5. We now present an alternative characterization of Sq{v) 
that may appear more intuitive to some readers. When J2j&i — Oi claim 
that the set So{i') defined in (4.11) also admits the following alternative 
characterization: Sq{i') is the unique set S that satisfies 

(4.12) ^(a^._^.)>o 
and 

(4.13) 7j < Oj + ^ ^ ("fc ~ 'yk) if and only if j G S. 

Since the strictly subcritical classes receive no extra capacity from the other 
classes, a necessary condition for the set S to be strictly subcritical is that 
the sum of the nominal capacities available to all classes in S is strictly larger 
than the sum of the mean arrival rates of classes in S — this leads to the first 
condition (4.12). Moreover, it is precisely this excess that is redistributed to 
the remaining classes, with class j receiving a fraction proportional to 
for j ^ S. This leads naturally to the "only if" part of the condition (4.13), 
which simply states that for any class that is not strictly subcritical, the 
total (reallocated) service capacity available to it is no greater than its mean 
(long-run) arrival rate. The "if" part of (4.13) is not as straightforward to 
justify a priori. Nevertheless, as shown below, it turns out to be the correct 
additional condition that uniquely characterizes Sq{v). 

To see that (4.12) and (4.13) uniquely characterize Sq{v), let S be any set 
satisfying (4.12) and (4.13). Since J2j£i^j = Jljeiilj ~ (^j) ^ Oi the relation 
(4.12) ensures that S ^2 and thus the relation (4.13) is well defined. Now, 
the two conditions above are equivalent to the relations 

(4.14) rs = > 
and 

(4.15) ^<rs iSjeS. 

Pj 

We now show that uniqueness of the set S is essentially equivalent to unique- 
ness of the representation (4.5). Suppose that we are given a set S^T 
satisfying (4.14) and (4.15). Define 

~ej = -{l- f3j)uj + - Pj)rs for j e S. 
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Then elementary algebra shows that (4.14) and (4.15) imply that 9j > for 
j G S*, and that R = v + J2jeS^kdk satisfies R = it{i^). However, uniqueness 
of vr dictates that R = k. Uniqueness of the representation (4.5) then shows 
that 9 = 9, with as in (4.11), and hence that S = So{i'). 

Now, let the sources be ordered so that 
(4.16) ^>^>...>i^. 

It then follows immediately from Remark 4.5 — in particular, relation (4.15) — 
that when J2jeJ^j ^ either So{iy) = or 5o(i^) has the form {j^, ■ ■ ■ ,J} 
for some j* SX. Given u G R"^, we now define the set of strictly subcritical 
classes to be 



(4.17) S{u) 



I, if 5^i^j<0, 

•So{i^), if XI - 0' 
i6J 



and let A4{i') be defined as follows: 

(4.18) M{iy) = {x€R-l:xi = for every i£S{iy)}. 

Note that M{u) = {0} is 0-dimensional if J^jex^j < M{u) = M.-^ is 
J-dimensional if X^jeJ^j — ^^'^ miuj vj > 0. The next result shows that 
A^(z^) acts as an invariant manifold for the fluid limit. 

Theorem 4.6 (Invariant manifold for the fluid limit). Given v G M*^ and 
u G '^l^, let K = 7t{i'), U = T{u + ui), and let S{v) and Ai{iy) be defined as in 
o-f^d (4-18), respectively. Then the following properties are satisfied: 

(i) //n G M{v), then U{t) = u + Kte M{iy) for all t G (0, oo). 

(ii) Given any compact set G C "^i^, there exists T = T{G) < oo such 
that for every uGG, U{t) G M{i^) for allt^> T. 

(iii) There exists u G Al(z^) such that Ui{t) > for every t G (0, oo) and 
ieI\S{iy). 

(iv) When J2jei^j —0? admits the equivalent representation 

(4.19) M{u) = {ueRi: T{u + ui){t) = u for all t G [0, oo)}. 

(v) Furthermore, ifJ2j£i^j = ^'^^ minjgji/j < 0, then there exists G 
I such that 

(4.20) ^ = ^ forj<j, and 5(z.) = |j G T: ^ < ^ j. 
Pj Pi I Pj Pi) 

Also, 

(4.21) z^j = - — forj^S{u). 
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Proof. Fix v G M-^. When Y^i^x^i < 0> '5(z>) = J, M{u) = {0} and so 
the first two statements of the theorem foUow directly from Lemma 4.4 of 
[17] (also see the results of [3]), and the third statement holds trivially. For 
the case when X^iej^i — minigxfi > 0, the first two statements of the 
theorem are a trivial consequence of the fact that = and the third 

statement is satisfied by any u with -Uj > for all i ^I, since for such u, 
T(u + vl) = u + ui. 

Therefore, for the rest of the proof of properties (i)-(iii), we shall assume 
that J2iej'^i ^ and min/giM/ < 0, in which case 5(z/) = So{i'). We start by 
proving property (ii). Given any compact set G C M.'^, let u* G be such 
that for every z G T, = max^gG Ui and let U* = T{u* + ul) . The comparison 
principle in Theorem 3.2 then guarantees that for any uG G, Ui{t) <U ^{t) 
for every z G X, where U = T{u + ui). Thus, in order to establish (ii), it 
suffices to show that there exists a finite T < oo such that U j (t) = for 
every j G So{i') and t >T. By the monotonicity property of the GPS SM 
established in Lemma 4.4 of [17], we know that 

(4.22) {i£l:U*{t)=0} C{i£l:U*{t + s)=0} forans,t>0. 
As a result, there must exist £ CI and T < oo such that 

(4.23) £ = {j(zI:U*{t)= 0} for all t > T. 

Since U* is piecewise affine, this implies that the slope of U*{t) for all t > T 
is well defined and equal to x> where x satisfies Xj = for j G £" and Xj ^ 
for j ^ £. On the other hand, by uniqueness of the SM, it is easy to see that, 
for s > 0, 

(4.24) U*{T + s) = T{u* + ui){T + s) = T{U*{T) + ui){s) = U*{T) + xs. 
From the definition of the SP, the last equality implies that for every s>0, 

{x-u)e d{U* (T) + xs) = cone[di,i G f] C d{x). 

Since x G I^+i by uniqueness of the projection vr, this show that x = ^ = ^(^)- 
The definition (4.11) of 5o(z^) and relation (4.6) of Lemma 4.4 then shows 
that 5o(z^) Q which completes the proof of property (ii). 

The definition of So{i^) implies that k — v £ cone[di : i G 5o(i^)] ^ cone[di : 
= 0] . If n G (z^) , then this implies that k — vCL cone[(ij : = Uj = 0] , and 
thus, by (4.5) and (4.6) of Lemma 4.4, f7(t) = T{u + i'i){t) = u + nt for ah 
sufficiently small t (and, in fact, for all t since k G K+). The above argument 
also shows that if n G M{u) with Uj > for every j ^ So{i'), Uj{t) = uj + 
Kjt > for all i > and j ^ So{u). This establishes properties (i) and (iii). 

Next, note that property (ii) shows that the set of invariant points of 
the fluid limit [described precisely by the right-hand side of (4.19)] must 
be contained in Ai. If, in addition, J2j£i^j ^ Oi then k = 7r(i/) = 0. Thus, 
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property (i) shows that u £ M imphes U{t) = T{u + i'i){t) = u for every 
t G [0,00), which estabhshes the fourth property. In particular, this imphes 
that X = = 0. Relation (4.20) is then an immediate consequence of relation 
(4.8) and the fact that 1 ^ Sq{v) [since Sq{i') t^T]. In turn, (4.20) implies 
that Vj = (3jUi/l3i for j ^ Sq^v). Summing the last equality over j ^ So{i^), 
and using the fact that Sfcex^fe ~ ^ J2keiPk = 1, we obtain = 
- J2k(^Soiu) ^k/ (1 - EfcG5o(i^) Pk)- When combined, the last two relations yield 
(4.21), thus completing the proof of the theorem. □ 

Remark 4.7. Although we do not use this later in the paper, it is 
also possible to provide a purely dynamical systems characterization of the 
invariant manifold M.{v) defined in (4.18). Suppose, given v S M'', an at- 
tractor A{v) for the fluid limit is defined to be any cone in M::[^ that satisfies 

(a) u G A{i>) implies U{t) = T{u + i'i){t) G A{i') for every t G [0, 00), and 

(b) u ^ A{v) implies \\m.t-^ooU {t) G A{i'). Then it can be shown that the 
invariant manifold M{v) is the intersection of M:j^ with the affine hull of any 
attractor A{u) for the fluid limit. In the subcritical case it is easy to see 
from Theorem 4.6 that the invariant manifold M{u) is the unique attractor. 
Indeed, as shown in Theorem 4.6(iv), in this case the invariant manifold 
M{u) can equivalently be characterized as the collection of invariant points 
for the fluid limit, which is the "standard" definition of an invariant manifold 
given, for example, in [2] and [13]. However, while the standard definition 
is limited to the subcritical case, our definition is more general in that it 
also applies to the supercritical case, where it provides information on how 
trajectories escape to infinity. Indeed, in the supercritical case it is not hard 
to show that a set is an attractor if and only if it is a cone contained in 
M{i') that contains the ray {nt^t > 0} and has nonempty interior relative 
to A^(z^). Since we do not use this property later, we omit a rigorous proof 
of this statement. 

5. Diffusion approximations for the unbalanced GPS model. For sim- 
plicity, we assume throughout this section that the classes are numbered 
so as to satisfy the ordering (4.16). As in the previous section, we consider 
a sequence of networks with associated processes H^,n G N, that satisfy 
Assumptions 2.1 and 4.1. Recall the defining equations for the fluid limit 
processes C/, X, Y and T given in (4.3) and (4.4), and consider the associ- 
ated diffusion scaled processes defined by 



(5.1) 



f/" = ^[C/"-f/], 
y" = y^[F"-F] 
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To prove the heavy traffic hmit theorem for the unfinished work process, 
we first assume that the primitive sequence {H^} satisfies, in addition, a 
functional central limit theorem. Let 7 € M.'^ be the vector in Assumption 
4.1(3) and let = 7 — a. 

Assumption 5.1. 

1. There exists a random variable u £ A^(i/) such that P a.s., 

lim — ■;=— = u. 

2. As n ^ 00, 

where B \s a, driftless J-dimensional Brownian motion with covariance 
matrix . 

3. There exists c G M-^ such that 

lim y/ni-^'^ — 7) = c. 

n — ^00 

We will also assume the heavy traffic condition 

J 

(5.2) E7^ = l• 

i=l 

Remark 5.2. (a) Note that Assumption 5.1(3) implies, in particular, 
that 7" — > 7 as n — > 00. 

(b) By the Skorokhod representation theorem (see, e.g.. Theorem 1.8 
of [10]), there exists a probability space {i},F,P), on which are defined 

P([0, 00) : ]R'^)-valued random variables H^, n G N, and B such that = 

H^^ B = B and H'^ B a.s. uniformly on compact sets (u.o.c). By an 
abuse of notation, we simply take Assumption 5.2(2) to mean that i/" B 
a.s. u.o.c. Of course, what we ultimately prove is weak convergence and not 
a.s. convergence, and this is reflected in the theorem statements. 

In the next section we summarize our approach to the GPS diffusion 
limit, and discuss its connection with related work. This section can be 
safely skipped without loss of continuity. 

5.1. General approach to the unbalanced GPS diffusion limit. The fluid 
limit result summarized in Theorem 4.2 shows that under the heavy traffic 
condition (5.2), [7 = 0. By (5.1), (2.5) and basic homogeneity properties of 
the SP, it then follows that 

[/" = ^/^" = r(V^) = r(V^), 
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where T is the GPS ESM and the last equahty fohows by Remark 2.4. If 
y/nX^ could be shown to converge to a limit X, then, since the GPS ESM is 
Lipschitz continuous, the continuous mapping theorem would immediately 
yield convergence of ?7" to T{X). Combining (2.7) with Assumption 5.1, it is 
not hard to see that \fnX converges if and only if 7 = a. This explains why 
the standard continuous mapping approach works only in the balanced case 
(see [17]). Indeed, in the (truly) unbalanced case, y/nX^ certainly diverges 
because the long-run average arrival rate 7^ for at least one critical class 
must be strictly greater than its nominal capacity. This class becomes critical 
(in the fluid limit) only because it receives extra service capacity from the 
strictly subcritical classes. Thus, the methods of [17], which carried out an 
analysis of the balanced GPS model, are not sufficient to analyze the more 
general, unbalanced case considered here. 

A similar generalization was considered in the context of single-class 
queueing networks. Specifically, the diffusion analysis in [19] [where each 
station was assumed to be in heavy traffic — see condition (24) therein] was 
extended in [4] to the unbalanced case. As mentioned in Section 4.2, the un- 
finished work U* for the model in [4] is represented, as in this paper, in the 
form U* = X* + y*, where X* is the associated netput process, but where 
Y* now admits the decomposition Y* = D^* for a nonsingular matrix D and 
1^* is the component- wise nondecreasing process that characterizes the cumu- 
lative idle time or excess capacity at each station. The ith column of D in [4] 
is analogous to the direction of constraint di in this paper, and the mapping 
from X* to U* corresponds to the SM T* considered here. The diffusion anal- 
ysis in [4] strongly uses (i) the explicit representation Y* = D^* = J2i=i Q'^i 
(see the block decompositions in the statement of Theorem 6.1 of [4]), as well 
as (ii) the continuity of the mapping that takes X* to which is referred 
to there as the regulator mapping (see (3.3D)~(3.3E) and the discussions 
following (4.25) and (4.29) in [4]). 

On the other hand, no such explicit representation is available for the GPS 
SP — in fact, the directions of constraint are linearly dependent, as shown 
in (3.1). Nevertheless, as shown in Lemma 4.4, an analogous decomposition 
holds when the GPS SM acts on affine trajectories. Indeed, suppose V is 
an affine trajectory, say, of the form u + vi for some u € M;|^ and v € M*^, 
and (j) = r(^), where F is the GPS SM. Since (j) is piecewise affine (i.e., 
with a finite number of changes of slope), from Lemma 4.4 and the property 
(4.24) of the SP one can infer that there exist unique measurable functions 
6 : [0, 00) — > with the property that for a.e. s G [0, 00), 

J+i 

^{s)=i;{s) + Y,ej{s)dj, 
i=i 

where, for j G T, 

(5.3) {i:0,(s)>O}C{i:</,,(s) = O} 
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and 9j^i{s) > implies = 0. Another application of the property (4.24) 
shows in fact that this also holds for piecewise affine i/j. Let ^(t) = /q 9(s) ds 
for t S [0, oo), and let G be the mapping defined on piecewise affine, contin- 
uous functions that takes ip to the corresponding ^. In view of the methods 
used in [4] discussed above, it is natural to ask whether B can be extended 
to define a continuous mapping on the space of all continuous functions. The 
following lemma answers this question in the negative — in fact, it shows that 
there is no canonical way to extend the definition of G to all continuous tra- 
jectories. This constitutes a structural difference between the GPS SP and 
other SPs considered in the literature in the context of queueing networks. 

Lemma 5.3. The mapping cannot he extended to define a continuous 
mapping on the space of all continuous functions. 

Proof. Let r be the 2-dimensional GPS SM, with directions of con- 
straint di = (1,-1), (i2 = (—1,1) and d^ = (l,l)/\/2- If there exists a con- 
tinuous extension of G to all continuous trajectories, then by the Cauchy 
property, for any pair of sequences of piecewise affine continuous functions 
{V'^^'^^lngN and {-(/-(^'"^InGN such that - V^^'"^ ^ 0, we must have 

0(^^(1,"-)^ — G(^/'^^'"'^) 0. We show that this is not the case by constructing 
a counterexample. Define ^(^'") = 0, let ^'-^'"^(0) = 0, and let 




if t G [2{m - 1)2-", (2m - 1)2""), 
if t G [(2m- l)2-",(2m)2-'^) 



for m = 1, . . . , 2" . Then it is easy to see that 

sup |^(i'")(s)-V(2,n)(^)|^ |V'(2'")(s)| < V2/n, 
se[o,i] «G[o,i] 

and thus, V'^^'") - V'^^'"^ ^ 0. On the other hand, if ^(i'") = 0(V'(^'")) and 
^(2,n) ^ 0(^{2,n))^ g^^g^jy verified that = 0, while ^^^'"^(l) = T^'^/n 

for i = 1,2 and C3^'"^(l) = 0. This implies that 

lim |C(2,n)/iN _^{l,n)(jN, 
n — ^oo 

which completes the proof. □ 

The aim of the above discussion was to explain why the analysis of the 
unbalanced GPS model does not fall into any of the previously existing 
frameworks for establishing diffusion approximations of queueing networks. 
We now briefly describe the approach taken in this paper to establishing dif- 
fusion limits, which entails first showing that there is state-space collapse. 
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in the sense that the strictly sub-critical classes vanish in the diffusion limit, 
and then using this information to provide a nice characterization of the be- 
havior of the remaining, critical classes. Lemma 5.3 suggests that it would be 
preferable to work, as far as possible, in the pre- limit since the GPS mapping 
is better behaved on the pre-limit than on the limit functions. Thus, we first 
introduce a modified work arrival process in Section 5.2 and identify the dif- 
fusion limit of the unfinished work associated with this work arrival process. 
The modified work arrival process is second-order equivalent to the original 
one in the sense that their diffusion limits coincide (see Theorem 5.4). How- 
ever, as shown in Lemma 5.5, the advantage of working with the modified 
arrival process is that the state-space collapse takes place in the pre-limit 
itself. This facilitates a simple characterization of the critical classes in terms 
of a certain reduced map, which is introduced in Section 5.3. The reduced 
map can be characterized as a GPS SM on a lower-dimensional space (as- 
sociated with the critical classes) with appropriately modified weights and 
is thus continuous. These results are finally combined with the comparison 
principle of Section 3.2 in order to establish state-space collapse and obtain 
an explicit characterization of the diffusion limit in Section 5.4. The diffu- 
sion limit identifies precisely how the covariance structure of the unfinished 
work of the critical classes is influenced by the variance of the cumulative 
work arrival processes of the strictly subcritical classes. 

5.2. A modified work arrival process. In this section we introduce a se- 
quence of modified cumulative work arrival processes {7i"',n £ N} obtained 
by smoothing the sequence of class j arrival process {-ff", n G N} for j £ S{i') . 
In Theorem 5.4 below, we show that {7{^,n € N} is equal to {H^,n € N} 
up to second order in the sense that their diffusion limits coincide. This is 
convenient because, as shown in the next section, the limit of the sequence 
of unfinished work processes associated with the modified arrival processes 
can be obtained by applying Theorem 5.6. 

A rigorous definition of the sequence {Ti^^n G N} is given below. Let 
{e",n € N} be a sequence of positive numbers such that 

(5.4) hm = and lim ^/ne^ = oo. 



+00 n^oo 



For n € N, define 

;,n^f7; + e", forje5(z.), 
~\7", for 3£l\S{u) 

and = 7" - a. For j ^ S{u), let H] = and for j G S{v), let H] be the 
departure process from a queue that is initially empty, has cumulative work 
arrival process H'^ and a deterministic service rate of 'yj'. In other words, 
we can write 



(5.6) = I 



/7;-ri(/7;-7^), foriG5(z.) 



^7, fori^5(z.) 
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where Fi is the one-dimensional SM defined in (3.9). It should be noted that 
TC^ does not satisfy Assumption 2.1(2). Nonetheless, it plays an essential 
role in proving our heavy traffic limit theorem. Let the corresponding fiuid 
scaled and diffusion scaled processes and W be defined as in (4.1) and 
(5.1), respectively, with H replaced by Ti.. The following theorem states a 
second-order equivalence result between the sequences {H"^} and {TC^}. 

Theorem 5.4 (Properties of the modified arrival sequence). Suppose 
Assumptions 4.1 and 5.1 are satisfied. Then the sequence {TC""} satisfies the 
following properties: 

1. n]{t) - n]{s) < ^f{t - s) for anyO<s<t and every j G S^v). 
1. Ti^^ 7""/- u.o.c. as m ^ oo, where 7i'^^{t) = . 
3. =^ -B as n ^ oo. 

Proof. From (5.6) and the explicit expression (3.9) for the one-dimensional 
SM we see that for t G [0, oo) and j G 5(1^), 

(5.7) rq{t) = - sup [7^5 - H^s)] V 0, 

se[o,t] 

from which the first property immediately follows. Next, note that by (5.6), 
the definitions (5.5) of 7 and (5.1) of ff", and elementary scaling properties 
of the one-dimensional SM, we have 

Tff = Hf - Pi (Hf - 7^ - e^'O and = - Pi {H^ - y^e^O • 

Taking limits as m ^ 00 in the first equation above and using the first limit 
in (5.4), along with Assumption 4.1(2), the fact that Pi is continuous and 
Pi(0) = 0, we obtain the second property of the theorem. 

By Peterson's "crushing lemma" ([15], Lemma 2), Pi(i/" — yJne^C) 0. 
[Although Peterson assumes that > e > 0, his proof goes through with 

— > if y^e" — > 00. Also, Peterson assumes mutually independent i.i.d. 
sequences of inter arrival and service times, but his proof goes through un- 
der our Assumption 5.1(2).] The third property is then an immediate con- 
sequence of Assumption 5.1(2) and Remark 5.2. □ 

For n G N, define 

0, ifje5(z/). 



^;(0)-|f/«(0), if j^S{u) 

and let A'" = Z^"(0) W - at. Moreover, let Z^" = F(A'"). 

The next lemma identifies the property that makes the unfinished work 
associated with the modified arrival sequence {TC^,n G N} easier to analyze 
than that associated with the original arrival sequence {H"',n G N}. Indeed, 
it establishes state-space collapse for the pre-limit sequence of unfinished 
work processes associated with the sequence of modified arrival processes. 
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Lemma 5.5. Suppose Assumptions 4.1, 5.1 and the overall heavy traffic 
condition (5.2) are satisfied. Then for every uj £^1, there exists N = N{u) < 
oo such that, for every n > N{uj), U'^{t){uj) G Ai{u) for all t G [0, oo). 

Proof. Fix w G and drop the explicit dependence on uj. Recall that 
= 7" - a, and define = r(Z^"(0) + and = r(Z^"(0) + vi). Then 
Theorem 5.4(1) along with the comparison principle of Theorem 3.2 shows 
that U^{t) < l>'"(t) for t G [0,oo). Thus, to prove the lemma, it suffices to 
demonstrate the existence of < oo such that, for n> N, <I>"(t) G 
for every t G [0, oo). By Lemma 4.4, there exist 5" > and x"' ^ such 
that, for f G [0,<5"), 

|."(t)=^"(0)+x"t and r'(t)=Z^"(0)+x"t. 

In addition, for every n G N, since U"'{0) G Theorem 4.6 shows that 

x" = K = 7r(i^). Moreover, for n G N, let 6'" G M."^^ be the unique vector 
in the representation (4.5) for x" and let G M:^''"^ be the corresponding 
unique vector in the representation (4.5) for k. Lipschitz continuity of the 
GPS SM r and the convergence — > for n — > oo, which is guaranteed by 
Assumption 4.1(3) and the first relation in (5.4), imply ^ k as n — > oo. 
The uniqueness of the representation (4.5) then dictates that On ^ as 
n — > oo. Thus, there exists N <oo such that, for every n> N, 

Soii^) = {j G T: 0,- > 0} C {j Gl-.e^ > 0} C {j G J:^;(0) = = 0}, 

where the first equality follows from definition (4.11) and the last inclusion 
is a result of relation (4.6) of Lemma 4.4. This implies that Xj = for 
j G So{i^) = S{i'), so that ^"(i) G Miiy) for t G [0,5"'). By the monotonicity 
property of the GPS SM, which was proved in Lemma 4.4 of [17] and is 
spelled out in (4.22) here, it then follows that $j(t) G Al(z^) for every t G 
[0, oo), which completes the proof of the theorem. □ 

5.3. A reduced representation for the SM. We assume throughout this 
section that 5 is a subset of I, S ^ I. The application of the main result. 
Theorem 5.6, of this section will be to the case when S = 5(z^) is the set 
of strictly subcritical classes associated with some such that X!j=i = 0- 
However, the results of this section are valid for an arbitrary set S C I. 
In order to formulate our reduced representation for the SM, we need the 
following notation. Let K, = X\ S and K = |/C| > 1. Given a weight vector 
[5 G (0, 1)^ (with EiLift = 1), let (3^ = {I3^,k£ K) be the JT-dimensional 
vector defined by 
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Also, for i ^ }C, let ef be the i('-dimensional unit vector associated with 
the ith coordinate, which has 1 at the ith coordinate and for all other 
coordinates in /C. Since /C / and /? G (0,1)"^, it follows that J^keicf^k > 
0. Hence, is weh defined, /3'^ G (0,1)^ and Y^ksK^k = which shows 
that is a K- dimensional weight vector. We now define the GPS SP 
associated with the weight vector f3^. First, let the directions of constraint 
{d^^,i G /C,c?^_|_i} be defined in terms of the i^-dimensional weight vector 
in the same way that the directions of constraint {di,i gX, dj+i} were 
defined in terms of the J-dimensional weight vector P in (2.3): more precisely, 
let 4^, = ^.g^ef/\/K and 

(5.9) df=ef- E foriG/C. 

Recalling definition (5.8), for i G /C, we can then write 

f 1, if j = h 

(5.10) {df), = , 

I EfcG/C Pk- 

Define 



if iG/C\{i}. 



A^ = {{xi,i G /C) : G M for ? G /C}, 

A^ = {{xi,i£lC):xi>0 for ielC} 
and, analogously to the definition given in (2.4), let 

d'^(x) = \ V Uidf : ffli > for i G l'^{x) 




where, for x G A^, l'^{x) is defined as in (2.5), but with T replaced by JC and 
J by K. Similarly, let the "/C-reduced" GPS SP be as described in Definition 
2.2, but with M-^, M;^ and d{-) replaced by A^, A^ and d'^{-), respectively. 
Finally, let F'^ be the associated SM, and note that it is well defined for 
the same reasons given in Remark 2.4. We now establish a correspondence 
between the original GPS SM and the /C-reduced SM. 

Theorem 5.6. Given a weight vector (3 G (0, 1)"^, let F he the associated 
GPS ESM, kt'0GP([O,oo);M'^), let = T{'il}) and for t e [0,oo) , letE{t) = 
{i ^I:fi{t) = 0}. Moreover, suppose there exist < Tq < T < oo and S C^X 
such that S C E{t) for every t G [Tq^T). Then, if K, = Z\S , the vector-valued 
process ip^ = {ipi,i G /C} satisfies 

ip'^{t)=T'^i^^)it) fortG[To,T), 
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where T is the fC-reduced GPS ESM defined above and ip'^ = £ /C} 

is given by 

(5.11) = V>i + /3f (^V^,) forieJC. 

Finally, if e dom(r), then ip^ = T^{i;'^). 

Proof. Since 1C^0,P^ and r are well defined. The theorem follows 
trivially if /C = T. Thus, throughout the proof, we assume that /C C T or, 
equivalently, that 5/0. Since T is single- valued on its domain (see Re- 
mark 2.4), if 99 = T{ip)^ then (p = r('0) for t € [0, 00), where 

^{t) = ^{TQ + t) and ^(t) = ^(ro + t)-v(ro) + y^(ro). 

Thus, we can also assume without loss of generality that Tq = 0. 

For every t £ [0, T), <f{t) € and, thus, (p'^{t) £ A+. Define r] = (p-ip and 
rj^ = — Tp^, It remains to show that in the interval [0, T), {(p^,r]^) satisfy 
properties 3 and 4 of Definition 2.3 for the GPS ESP on (corresponding 
to the weights (3^). First, note that for i G /C and t € [0, T), 

(5.12) 

= 7?,(t)+/3f (^%(i)), 
\k&S I 

where the second equality follows from the fact that ^p^^ (t) = for k £ S and 
(p = ilj + r]. 

Now, fix < s < t < T. In analogy with the hyperplane H defined in 
Lemma 3.1, let = {x G A'^ : J2ieic Xi = 0}- Recah that for every x € \ 
{0}, d{x) C H C d{0) = {x G M"^ : Y.i=i > 0} and, likewise, for every x £ 
A^ \ {O'^}, d'^ix) (ZH^C d^(O^) = {x £ : Y^ieic^i > 0}- Also, note that 
for any u £ [0,T], J2ieic V'f (u) = Y.iex^i{u) since ipi{u) = for i G 5, and 
also that (5.11) implies J^ieK.'^i^i'^) — i^ii'^) ■ The last two statements 
along with the fact that {<p, rf) satisfy property 3 shows that 

(5.13) Y.i^f{t) - r,f{s)) = Y,{r,,{t) - ^,{s)) > 0, 

with the inequality being replaced by an equality if, for every u £ {s,t], 
tp{u) / 0. This shows, in particular, that r]'^{t) — rf^{s) £ d^{Qi^) and rf"{t) — 
rf^{s) £ if, for every u £ {s,t], (p{u) / 0. Thus, this proves property 3 in 
the case when either ip{u) = for some u £ {s,t] or when ip{u) 7^ for every 
u £ is,t], but co[U„e(.,t] d'^ip'^iu))] = H'^. 
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We now consider the remaining case when <p{u) 7^ for every uG {s,t], and 
^[U«G(s t] d'^{ip^{u))] C H^. By Lemma 3.1, the latter relation implies that 
there exists j S /C such that, for every u S {s,t], ipj{u) = ^f{u) > 0. Since 
(if, rj) satisfy property 3 of the ESM and the vectors di,i £ E = I\ {j}, are 
linearly independent, there exist unique 9k>0, k £ E, such that 

(5.14) rj{t)-r]{s) = Y.ekdk 

k&E 

and 

(5.15) 0k>^ =^ V'kiu)=0 for some li G (s, t]. 

Substituting r/(t) — 77(5) for w in (3.2)-(3.3), and using the fact that 5 C 
the relation (5.14) implies that, for i£2, 



(5.16) ^ 
Summing over i £ S, this yields 

Uvdt) - mis)) = (e f 1 (e a) ( E 



The last two relations can be combined with (5.12) to show that, for i G /C, 



^f(()-^f(.) = -r\-ft(l + #^ 

1 - Pi V 2^keK Pk ^ 



1 - A " 

Using the definition 



(5.17) = ^ ^^ Qi forieJC 

1 Pi 



and the fact that E\S = E nJC, the last equation can be rewritten as 

for i e /C. 



vHt)-^hs) = -^-P^( E 
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From the definition of the /C-reduced directions of constraint given in (5.9), 
arguments analogous to those used to obtain (5.14) and the first hne of 
(5.16) can be used to show that 

keEnic 

That property 3 of the ESP is satisfied can now be inferred from the fact 
that, for k€lC, (5.17) shows that 6*^ > if and only if 6*^ > 0, and (5.15) 
ensures that dk > implies (pk{u) = [equivalently, 9?^(tf) = 0] for some 
u e {s,t\. 

The fourth property of the ESP can be proved in an analogous fashion, 
splitting into the cases when ip{t) = and ip{t) ^ [in which case B. ^ d{ip{t)) 
and the set {di : ipi{t) = 0} is linearly independent]. The argument is exactly 
analogous to the proof of property 3 and is thus omitted. Thus, we have 
shown that = r'^('0^). 

li ip ^ dom(r), then, by Theorem 1.3(2) of [16], is of bounded varia- 
tion on every compact time interval, which in turn implies that 77^ is of 
bounded variation on every compact time interval. Since {ip'^,r]^) satisfy 
the /C-reduced ESP for once again invoking Theorem 1.3(2) of [16], this 
shows that {(p^ ,7]^) in fact solve the /C-reduced SP for -0'^, thus establishing 
the last statement of the theorem. □ 

5.4. Heavy traffic approximation for the unfinished work. Fix v E R'^ 
with X^jeJ^i — a-iid denote the set of strictly subcritical classes S{v) sim- 
ply by S. Recall that, due to the ordering (4.16), S has the form {j*, . . . , J} 
for some j* > 1 (with j^, set to J -|- 1 if 5 is empty) . Let K = — 1, so that 
K is the cardinality of the set /C = X \ 5 of critical classes. It will be con- 
venient to introduce the linear "projection" operators £ : ^([0, 00) :M'^) — > 
2?([0, 00) : R^^) and L : R-^ ^ defined by 

[cfh = h + PnT.fi] if^e/c, 

VjG5 / 

for / G ^^([O, 00) :R'^) and, analogously, 

[Lv]^ = v^ + pHJ2vj] if^e-^ 

VjG5 / 

for V € R"^. (In all cases, the sum over an empty set is taken to be equal to 
zero.) Define = jO{u + B + cl), and note that (since u G A4) 

= Ui + Bi + CiL + P^T (^i + • 
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Moreover, recalling the definition of the /C-reduced GPS ESM T associated 
with the weight vector defined in (5.8), let 

and define the J-dimensional process U = {U^,0^), where here 0"^ is the 
identically zero process in the space A;^ = {{xj^ , xj) : x G M"^}. In other 
words, let 



(5.18) 




for i G /C, 
for i £ S. 



Also, define 

f = u + B-U + CL. 
We now state and prove the main result of the paper. 



Theorem 5.7. Suppose Assumptions 2.1, 4.1, 5.1 and the overall heavy 
traffic condition (5.2) are satisfied. Then ?7" ^ U and ^ T. 

Proof. For n G N, define W^'*" = A'^" = CX"", ZY^-" = r'^(A'^'^), 

B'^ = CB and, likewise, = La, = Lu, = Lc, = Lu and vF = Lu. 
Note that, since Assumption 4.1(1) and Assumption 5.1(1) together imply 
that u = and n G TW, it follows that u'^ = and uf = Ui for i G /C. By (4.20) 
and (4.21), we know that, for i G /C, 7,, - ai = Vi = -PiY.jes'^j/T.keicf^k- 
By (5.2), J2j=i — 0- When rearranged, this shows that = or, equiv- 

alently, that 7^^ = a^. Also note that J2i=i Q^i = 1- Now let U ' and X ' 
be the fluid scaled versions of W'^'" and X^'"^, as defined in (4.1), and also 
define = 0, = 0. Then, by Assumption 4.1(1) and the definition of 
Uj'{0) given after the proof of Theorem 5.4, we have 



lim u'^'''(0)=l( limlTm) <l( \im ^^^] = Lu = u'^ =0. 

n— >oo \ n-— »oo / \ n-^oo ji ' 



When combined with Theorem 5.4(2), Assumption 4.1(3), the linearity of the 
operators C and L and the fact that = 0, this shows that x^'^ X^ = 
u.o.c. Since is the SM associated with the GPS model with basic and 
redistribution vectors and (3^, Remark 2.4 ensures that it is Lipschitz 
continuous. Since X and X are of bounded variation and thus lie in the 
domain of the SM, this allows us to conclude that Z// — > T^{X ) = = U 
u.o.c. 

Now define H'^'" and Z?'^'" as in (5.1), with U,X and H replaced 

by Z//^, X'^ and Ti^, respectively. The fact that X^ = l/(^ = 0, the scaling 
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properties of the SM and the fact that the ESM coincides with the SM on 
the domain of the SP (see Remark 2.4) show that zJ'C'" = r^(i''C'") and 

\ n 



Along with Theorem 5.4(3), Assumption 5.1(1) and Assumption 5.1(3), this 
shows that ^ u.o.c. The fact that = r'^(^'^'") and the Lip- 



schitz continuity of the GPS ESM T (which follows from Remark 2.4) 
then shows that Z^'^'" ^ u.o.c. Now fix w G ri. Since Lemma 5.5, along 
with Theorem 5.6, shows that there exists N = N[uj) < oo such that, for all 

n > N{uj), 

an, N fO' if I g5, 

- \ [r'C(;e'C."(L^))]i =Z^f'"(a;), if i G/C, 

we conclude that W^iuS) U (lo) for every a; G il. 

Furthermore, since U = U = and F is Lipschitz continuous, there exists 
a constant C < oo such that, for any T < oo, 

sup |C/"(t) -Z?"(t)| 

tG[0,T] 

= sup |I7"(t) -Z7"(t)| 

tG[0,T] 

= ^/^ sup |F(X")(t) -F(:r')(t)| 

*e[o,T] 

<C^/^ sup |X"(t) -:Y"(t)i 

te[o,T] 



C/"(0) ^/"(O) 



n vn 



+ sup \7r{t)-7{'{t)\ 

te[o,T] 



< JC7 max ^-=^ + C sup |i7"(t) - W"(t)|. 

iG[0,T] 

Taking limits as n — > oo. Assumption 5.1(1) and the fact that u£ M show 
that the first term in the last line of the display tends to zero, while Assump- 
tion 5.1(2) and Theorem 5.4(3) show that the second term tends to zero. 
When combined with the fact that U"" U u.o.c, this shows that — > U 
u.o.c, as desired. 

Finally, using equations (2.1) and (2.2), for n G N and t G [0, oo), we can 
write 



Tr{t) = uno)+Hnt)-ur{t). 
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Using the relation T = ^l, which foUows from Theorem 3.2, this imphes that 

T^it) = ^ +Tr:{t) -Tr,{t) + v^(7r - l^)t. 

By Assumption 5.1 and the diffusion hmit for the unfinished work just proved 
above, we then have 

(5.19) f"" ^f = u + B-U + ci, 

u.o.c. By Remark 5.2, the pathwise u.o.c. hmits must be replaced by weak 
convergence, and the proof of the theorem is complete. □ 

Remark 5.8. It is also possible to introduce additional primitives to 
describe the queue length Q", sojourn time and waiting time VF" pro- 
cesses associated with this model. Indeed, this was done in [17] where, un- 
der general assumptions on the first-order properties of these primitives (see 
Assumption 4.5 and Condition 4.8 of [17]), fluid limits for these processes 
were established (see Theorems 4.7 and 4.11 of [17]). In particular, it was 
also shown in Lemma 4.13 of [17] that when the heavy traffic condition 
SiLi 7i = 1 holds (and the initial conditions converge in a suitable manner 
to 0), the fluid limits of all these processes are identically zero. In addition, 
under general functional central limit type assumptions (see Assumptions 
4.12 and 4.16 of [17]), diffusion limits for these processes were deduced from 
the corresponding diffusion limits for the unfinished work and busy time 
processes. Indeed, it was shown in Theorems 4.18 and 4.19 of [17] that the 
diffusion limit Q for the queue length process satisfies Qi = fiiUi for i G X, 
where l/^i corresponds to the long-run average service requirement of a 
class i customer, while the diffusion limits V and W for the sojourn and 
waiting time processes, respectively, satisfy Vi = Wi = Ui/ji for i £l. These 
diffusion results were obtained in [17] under the balanced heavy traffic con- 
dition 7 = a. However, it can be shown that these results continue to hold 
even in the unbalanced case considered in this paper. Indeed, they can be 
deduced from the results of Theorem 5.7 by following almost verbatim the 
proofs of Theorems 4.18 and 4.19 in [17], with the only modification that 
a be replaced by 7 in those proofs. As a consequence, we omit a detailed 
exposition of these results. 
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